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The zero range potential is constructed for a system of two particles interacting via the Coulomb 
potential. The singular part of the asymptote of the wave function at the origin which is caused by 
the common effect of the zero range potential singularity and of the Coulomb potential is explicitly 
calculated by using the Lippmann-Schwinger type integral equation. The singular pseudo potential is 
constructed from the requirement that it enforces the solution to the Coulomb Schrodinger equation 
to possess the calculated asymptotic behavior at the origin. This pseudo potential is then used for 
constructing a model of the imaginary absorbing potential for the positron electron system. This 
potential allows to treat the annihilation process in positron electron collisions on the basis of the 
non relativistic Schrodinger equation. The functional form of the pseudo potential constructed in 
this paper is analogous to the well known Fermi-Breit-Huang pseudo potential. The generalization 
of the optical theorem on the case of the imaginary absorbing potential in presence of the Coulomb 
force is given in terms of the partial wave series. 
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I. INTRODUCTION 

The Dirac-delta potential was introduced by Fermi 
to describe the interaction of neutrons with atomic nu- 
cleus [J. Already at that time it was realized that only 
perturbative treatment of such a potential is possible [2[ 
due to the strong singularity of the delta function in three 
dimensional. 

One of the possible way to go beyond the perturba- 
tive treatment is the use of the concept of zero range 
potentials 2-4]. In this approach we impose the singular 
boundary condition 



(r,fc)~- 

r 



(1) 



as r — > 0. Here 
equation 



(r, k) is the solution of the Schrodinger 



[-A r + V(r) - k 2 ] 4>{r,k) = 0, 



(2) 



where A stands for the Laplacian, bold letters are used 
for vectors in R 3 like r and k, and r = \r\ and k = \k\. 
The potential V(r) here should be of a short range type, 
i.e. it should be less singular than 1/r 2 at the origin 
and should vanish faster than the Coulomb potential 1 jr 
as r — > oo [5]. The parameters a s and f3 s are not in- 
dependent and should be correlated in such a way that 
the wave function (j>(r, k) is the solution of Eq. © when 
r — > 0. As a result, the only independent parameter a s 
of the model can be chosen as a s = —a s /f3 s . In the case 
of V(r) — this parameter is the scattering length. 
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To enforce the boundary condition in an alternative 
approach [6, 7] an additional term is introduced into the 
Schrdingier equation 

[-A r + V(r) - k 2 ] </>(r, k) + XS(r)l3 s = 0, (3) 

where A = Aira s . Now j3 s can be obtained as the limit 



lim —rd)ir, k) 
r^o dr 



(4) 



This allows us to rewrite (|3]) as 7] 

[-Ar + V(r) + XW s (r) - k 2 ] <p(r, k) = 0, (5) 

where 



W s (r) = S{r)—r 
ar 



(6) 



,k) 



is a kind of "pseudo" potential. We note that all 
above expressions can be derived [5[ from the Lippmann- 
Schwinger equation 

o (r,fc)-A J dr' G(r,r',k 2 +iO)W s {r')i>(r' ,k), 

(7) 

where 0o( r , k) is the regular solution to the equation (O 
with A = and G(r, r', k 2 + iO) is the Green's function 
for the Hamiltonian — A r + V(r). The explicit form of W 
given in Eq. ([B]) is valid only for short-range potentials 
and it should be modified for long-range potentials like 
the Coulomb potential Q. 

Zero range potentials are the mathematically correct 
tools for describing contact interactions. However, in ap- 
plications the delta function is often used as a potential of 
contact interaction. A recent example of such an applica- 
tion of the delta function potential is its use in describ- 
ing the annihilation of positrons in atomic systems 
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The loss of particles in positron-hydrogen scattering due 
to the e + e~ annihilation can be simulated by an imag- 
inary absorption potential in the positron-electron sub- 
system. Arguments from perturbative QED Fu| suggest 
a three-dimensional delta- function potential [9j. In a se- 
ries of calculations of positron hydrogen scattering [9llTiT- 
[l3j this delta functional potential was smoothed out and 
the delta function was replaced by Gaussian with finite 
width. Thus the concept of zero-range potentials was not 
used. Nevertheless, a mathematically sound formulation 
of how the three dimensional zero range potentials in the 
Coulombic systems can be used is of definite demand. A 
way for such a formulation was outlined in Ref . [l4| • 

In this paper we continue developing the zero range po- 
tential concept in the presence of Coulomb interaction. 
We give here the detailed derivation of the pseudo po- 
tential for the two body Coulomb Hamiltonian. We show 
that, as in the short range case, the Schrodinger equation 
for this potential can be solved analytically with the help 
of a Lippmann-Schwinger-type integral equation. We de- 
rive from this solution the exact asymptote of the wave 
function at the origin and then the corresponding pseudo 
potential is constructed. We show that, if this potential 
is treated perturbatively, it can be taken in the simple 
form of the delta function as it was used in Q for the 
positron electron case. Finally, we apply the constructed 
zero range potential to the description of the positron 
electron annihilation by introducing the imaginary ab- 
sorbing pseudo potential. With this potential we solve 
the bound state problem deriving the positronium spec- 
trum as well as the scattering problem calculating the 
positron electron wave function and the scattering am- 
plitude. The annihilation cross section corresponds to 
the loss of the flux due to the imaginary absorbing po- 
tential. It is calculated by using a generalized optical 
theorem which we have modified in such a way it can be 
applicable to the case of Coulomb potential. 



II. PSEUDO POTENTIAL FOR COULOMB 
HAMILTONIAN 



In this section we construct the zero range potential 
and the corresponding pseudo potential for the Coulomb 
Hamiltonian Hq 



where 



H c = -A r + n/r, 



2ne 2 Z 1 Z 2 
1 K 2 



(8) 



(9) 



equation 

(H c ~k 2 )^ c (r,k) = 0. 
The explicit form of ipc(i*,k) reads [l5| 



(10) 



ip c (r, k) = T(l + iT])e-^ 2 e ir - k M(-ir], 1; i[rk - r • fe]), 

(11) 

where r/ = n/(2fc), T(z) is the gamma function and 
M(a, 6; z) is the re gula r Kummer (or confluent hyper- 
geometric) function |16| . 

Now we construct the zero-range potential and the cor- 
responding pseudo potential along the scheme given by 
Eqs. (HJTJ). The modified Schrodinger equation reads 

[H c - k 2 ] i>(r,k) = -\W{r)^{r,k) (12) 

and W is such that 

W(r)i)(r,k) = S(r)/3. (13) 

The parameter /3 can be determined from the solution of 
the Eq. (TT2"j) . This solution, as in the short range case, can 
be obtained by solving the Lippmann-Schwinger equation 

ip(r,k) = <ip e (r,k)-\ / dr' G c (r,r',k 2 +iO)W(r , )tP(r',k). 



(14) 

Here Gq stands for the Coulomb Green's function. In 
view of Eq. (fT3"]) the integration on the right hand side of 
Eq. (fT4"|) immediately leads to the representation 

ip(r, k) = ip c (r, k) - \G c {r, 0, k 2 + i0)/3. (15) 

From the explicit form of the Coulomb Green's function 



given in 17[ one can derive the following representation 



in terms of the Whittakcr function 

.16] 

G c (r,0,fc 2 + i0) = -^ r ( 1 + i 'nW_ lri .i(-2ikr) . (16) 

This representation is the basis for computing the asymp- 
totic expansion of ip(r, k) as r — > 0. The Coulomb wave 
function is regular at r — and can be expanded as 

ip c {r,k) = T(l + ir 1 )e-^/ 2 [l + 0(r)}. (17) 

The Coulomb Green's function Gc(r,0,k 2 + iO) is sin- 
gular as r —¥ 0. We will determine the explicit form of 
this singularity by using the representation (|16l) . The 
Whittakcr function can be represented by the irregular 
Kummer function 11611 



-A* A 

e 2 2 



z^[/(l/2 + /i- 7 ,l + 2^,z), (18) 



which is valid for — tt < arg z < ir. With 7 = — ii] it 
results 



W 



(-2ikr) = -2ikre lkr U(l + ir), 2, -2ikr). (19) 



it]; 



with e = ±1 for the repulsive and attractive case, respec- 
tively. The Coulomb wave function V'c( r : k) satisfies the 



We use following representation of the irregular Kummer 
function expansion [16] for integer m 
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(_l)m+l 

U(a, m + 1, z) = — 7=-, r \M{a, m + 1; z) lnz 



E „/ a !'n *' W a + - ^(1 + - ^(1 + m + I)} 



m\T(a — m) 
(m — IV 

* — —-±z- m M{a -m,l-m; z\ 
T(a) 



(20) 



I 

Here — 7r < arg z < n, (a); = r(a + V)/T(a) is the The regular Kummer function M(a, b; z) is given in terms 
Pochhammer symbol, ip{z) = F'(z)/T(z) is the digamma of M{a, b; z) m by the limit 



function and M{a, b; z) m is the truncated regular Kum- 
mer function defined by the polynomial 



M(a, b; z) = lim M (a, b; z) 7 



M(a,b;z) m = ^2^z l 



1=0 



H(b)i 



(21) 



Using now (|2U1) for the U factor in (ITO1) with appropriate 
arguments and keeping the leading terms we arrive at the 
following asymptotic representation for that U factor as 

r ->• 



U(l + ir],2,-2ikr) 



1 



1 



2ikY{l+irf) 



[1/r + 2fc^log7 



17] 



nT(l + it}) 



[4irC(k) - ik] +C(rlogr). 



(22) 



Here C(k) is given by 

C ( k ) = ^ + ^ [ l0 S(- 2lfc ) + i> (! + + 2 7o - 1] , 

n (23) 

where 70 is the Euler-Mascheroni constant [161 ] . Sub- 
stituting (|19p and (|22l) into (jT6j) we obtain the required 
Coulomb Green's function as r — >• 

G c (r,0,fc 2 -MO) = — [1/r + nlogr] +C(fc) +Ofrlnr). 

47T 

(24) 

With these results introduced into (|T5|). we end up with 
the following asymptotic expansion for the wave function 
as r — > 



V>(r,fc) = 



-A/3 

47T 



[1/r + 11 log r] +Vreg + 0(rlogr), (25) 



where "0 rO g = T(l + irf)e^^' n / 2 — \f3C{k) is the regular 
part, which can be idetified as j3 — ip lcg - The asymptotic 
expansion of the wave function as r — ¥ now takes the 
final form 

ip(r,k) = — [1/r + nlogr] + /3 + 0(r logr), (26) 

4-7T 



and the constants are given by 
a — —A/3, 



T(l + i?7)e 



-7T7J/2 



1 + AC(fe) 



(27) 
(28) 



As in the short range case, the only independent param- 
eter A now determines the zero range potential. Similar 
to the short range case A is related to the ratio of param- 
eters a and j3 by A = —a/f3. 



I 

At the last stage of our construction let us represent /3 
in terms of the wave function. Introducing a new variable 



1 + nr log r 



it is easy to obtain by direct calculation the following 
result 



lim —uil)(r, k) = B. 
r^a du 



(29) 



This formula allows us to rewrite the pseudo potential W 
in the form 

W(r) = S(r)4-u , 
du 

and even more explicitly in the final form 
(1 + nr log r) 2 d r 



W(r) = S(r)- 



1 — nr dr 1 + nr log r 



(30) 



We can see that in the non-Coulomb case, when n = 0, 
we recover the short-range result of Eq. ([6]) . 



III. STRUCTURE OF THE GREEN'S 
FUNCTION 



In this section we construct the Green's function for 
the Hamiltonian Hq + \W(r). It is the solution of the 
equation 



[H c + XW — z] G(r, x, z) = S(r - x) 



(31) 
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The Green's function is symmetric with respect to r and 
x and obeys the boundary condition as r — > 

A(x z*) 

G(r,x,z) = V ' [l/r + nlogr]+B(x,z)+0(rlogr) . 

4.7T 

(32) 

Like in Eq. (ITBl , the action of the quasi potential on the 
Green's function is given by 

W(r)G(r, x, z) = 5{r)B{x, z). 

Thus, we can combine (|3"Tj) and (I3"2l to obtain 

[H c - z] G(r, x, z) = S(r - x) - XS(r)B(x, z). (33) 

By inverting the operator He — z with the help of the 
Coulomb Green's function, and performing the integra- 
tion with the the delta function, we arrive at the repre- 
sentation 

G(r,x,z) = G c (r,x,z) - XG e (r,0,z)B(x,z). (34) 

Taking the limit when r —¥ for a fixed nonzero x with 
the use of (|2~4"|) we arrive at the expression 

~, * XB(x 7 z) r . . 
G(r, x, z) = [1/r + nlogrj 

47T 

+G C (0, x, z) - XB(x, z)C(y/z) +0(r log r). (35) 



Comparing this asymptotic form with the definition 
we find the expressions for A(x, z) and B(x, z) 



A(x, z) — Gc(0, x, z)- 



-X 



1 + \C{y/z) 

B(x,z) = G c (0,x,z) 1 . (36) 
1 + XC(y/z) 

By substituting B(x,z) into (|34p . we obtain the final 
representation for the Green's function in the form 



G{r,x,z) = Ge(r,x,z) 

-XG c (r,0,z) 



1 + \C{y/z) 



G e (0,x,z).(37) 



As a straightforward consequence of this representation, 
the equation 



1 + XC(Vz) = 



(38) 



determines the poles of the Green's function, or in other 
words, the spectrum of eigenvalues of the Hamiltonian 
H c +XW(r). 

The Eq. (|37D can be written in the operator form 

G(z) = Gc(z)-G c (z)T(z)G c (z), (39) 

were the T matrix is defined by the kernel 

S(r)5(x) 



T{r,x,z) = X 



1 + AC(v^) ' 



(40) 



If now we formally introduce the Coulomb T-matrix by 
the representation of the Coulomb Green's function 

G c (z) = G (z) - Go(z)T c (z)G (z), (41) 

where G (z) = (-A-z)~ l , then the Eq. fl39]) can further 
be rewritten in the standard form 



G(z) = G (z) - G (z)T(z)G (z). 



(42) 



Here the T-matrix T(z) for Coulomb plus zero range po- 
tentials is given by the explicit expression 

T(z) = T c (z) + [I- T c (z)Go(z)]T(z)[I - G (z)T c (z)}. 

(43) 

This formula reproduces the conventional representation 
for the T-matrix in the case of the Coulomb plus short- 
range interaction which follows from the standard two 
potential formalism [l5j |. 



IV. ELECTRON POSITRON ANNIHILATION 
POTENTIAL 

The formula for the 27 singlet e + e~ annihilation cross 
section derived from perturbative QED reads EH 



farm = CoZ c ff, 



(44) 



where Co = irr^c/v), tq is the classical electron radius, c 
is the speed of light and v is the incident velocity of the 
positron. The effective number of electrons Z c g partici- 
pating in the annihilation is given by the integral 



Z cfi = J dr 1 dr 2 \^ 0+ {r 1 ,r2,r p )\ 2 5(r 1 -r 2 ) 



(45) 



Here \I/ 0+ (ri, r 2 , r p ) (with 7"i, r 2 and r p being the posi- 
tion vectors of the positron, the electron and the proton) 
is the solution of the positron hydrogen scattering prob- 
lem, when annihilation is not taken into account, and 
normalized in such a way that the incident wave has the 
form 



exp(ifci • ri)$(r 2 ,r p ), 



(46) 



where Q(r 2 ,r p ) is the hydrogen wave function. 

In the Ref. to describe the absorption in the 
positron hydrogen scattering process an imaginary op- 
tical potential of the form igW\ 2 , with g < 0, was in- 
troduced. The optical theorem, in the first order Born 
approximation, gives the following expression for the ab- 
sorption cross section 

B _(~g) 

ki 

x J d ri dr 2 W 12 ( ri - r 2 )|*°+(r 1 ,r 2 ,r p )| 2 .(47) 

By comparing the QED expression (j4"4")l with Eq. (|4T|) 
and taking into account that cr ann = o^ n /4 (the factor 



1/4 comes from the fact that only singlet state is to be 
taken into account) it was proposed in Q to determine 
the optical potential in such a way it has the coupling 
constant of the form 



g = -kia 
and the coordinate part of the form 
Wia(r) =S(r), 



(48) 



(49) 



where r = T\ — r 2 . 

As we have discussed above, such a choice of the coordi- 
nate part of the potential leads to troubles in applications 
since the delta function as a potential can only be treated 
on the basis of the perturbation theory. Moreover, as it 
was studied in details in Ref. [14j . the formula (|47p is 
not valid for energies above the positronium formation 
threshold. In that paper it was proposed to use the zero 
range potential as a remedy for all of these problems. 
Thus, the actual positron electron annihilation potential 
in the framework of above considerations should be de- 
fined as 



igW 12 (r) = igW(r), 



(50) 



where W(r) is the pseudo potential of Eq. (|30|) . Con- 
sequently, the positron electron Hamiltonian, which si- 
multaneously describes the inter particle dynamics and 
annihilation, takes the form 



H = H C + igW{r) 



(51) 



THE POSITRONIUM SPECTRUM 



In this section we study the positronium spectrum as 
the discrete spectrum of the positron-electron Hamilto- 
nian He +igW(r) with annihilation potential (|50"1) . Here 
is the case when e = — 1 in Eq. ^ and n < 0. The spec- 
trum of H is formed by zeros of the denominator in (|37[) 
when X — ig. For analyzing solutions of the equation 
we rewrite it in the form 



(52) 



where 



4tt 



\og{-2i^) + i> (\ + ^^=^ + 2 7o 



- 1 



(53) 



Since the actual value of the coupling constant g is small, 
around 10 -6 eV, we can seek the solution of Eq. (|52|) 
in the g — > limit. In this limit the right hand side 
of (|52p tends to infinity and so should do the left hand 
side. It may happen only due to the digamma function. 
Using the well known series expansion of the digamma 
function 03 



we see that if the digamma function argument is equal 
to — m, where to is a non- negative integer, the left hand 
side of Eq. ([52|) will obey the required property. Thus, 
we get the equation that determines the zero order ap- 
proximation to energy levels 



1 + i 



2^i 



TO = 0, 1, . 



(55) 



Solving it, we recognize the pure Coulomb spectrum 

T ,2 



ZN 



AT = 1,2,. 



(56) 



in accordance with the obvious expectation when no ab- 
sorption is present in the Hamiltonian if g = 0. 

As the next stage we will calculate the first order cor- 
rection to zjy. First we represent the exact energy level 
in the form 



z' N = -(xiv + <5jv) 2 , 



(57) 



where — x 2 N = zjq is the Coulomb energy and Sn is the 
correction of interest due to the presence of the annihila- 
tion potential. Obviously, we expect to have Sn — O(g). 
For calculating Sn let us use the inverse to the equa- 
tion ([521 



1 



-ig. 



(58) 



According to the discussion above, the digamma function 
argument for values of z in the vicinity of zjv tends to a 
nonpositive integer when g — > 0, i.e. 

l + i-^= = l-N--^6 N + 0(g 2 ) N = l,2,.... 



N 



(59) 

Due to Eq. (|54l) for the leading term of the digamma 
function when x — > —p with integer p we have 



ip{x) 



1 



x + p 

So, we can conclude that 
1 



P = 0,1,. 



(60) 



^ 1 + ^ 2 



2k 2 n 



S N +0(g 2 ), (61) 



where N — 1, 2, . . .. Now the left hand side of the equa- 
tion (l58l) takes the form 



1 



-in 



+ 0(g 2 ). (62) 



Considering (|61[) we can rewrite this in the form 



6 



and for Sn we obtain 

Sn : 



+ 0(g 2 ). (64) 



2vr ' ~ K3 ' 8ttN 2 
Consequently, the energy levels are given by 

(-n) 3 

z' N = -k 2 n - 2k n 5 n -S n = -k 2 n + ig 3 + 0(g 2 ) 

(65) 

with N = 1, 2, . . .. So we got the first order correction to 
energy levels due to the annihilation potential 



(66) 



It is interesting to compare this value with the result 
of the standard perturbation theory. For the potential 
igW(r) the first order energy level correction is given by 
the matrix element 

(ipNe m \igW\ip Nem ) = J dr ip Nem (r)igW(r)ip Nem (r), 

(67) 

where i^Nim(f) is the Coulomb bound state wave func- 
tion pi 



4>Ni m {r) = 



(-n) 3 2 (N-£-l)\ 



8 N 2 V [(N + £)\] 



r nr i / — nr 
XeXP l2lv} — ' 1 



21+1 
N-i-1 



N 



The integral on the right hand side of (|67|) can easily be 
evaluated to give 



(ipNhn\igW\ip N i m ) = igip 2 Nlm (0). 



(69) 



Since V'jvzrn(O) is non-singular only for s-wave, i.e. for 



Nl 

I = 0, we have 



^oo (0) 



(~n) 3 
8ttN 3 ' 



(70) 



This leads us to the final result for the perturbative cor- 
rection 



(ipNoo\igW\ip N Q ) = ig 



(71) 



which is in complete agreement with our previous result. 

VI. ELECTRON-POSITRON SCATTERING 

In this section we consider the positron electron scat- 
tering with annihilation. The latter is described by the 
annihilation potential igW(r). The corresponding wave 
function t/j(r, k) is the solution of the Schrddinger equa- 
tion (IT21 with A = ig and it possess the asymptotic be- 
havior as r — > oo 

ip(r, k) exp[ifc ■ r + i-q log(fcr — k ■ r)] 

| f^ exp[ikr - ir]log(2kr)] ^ 



where 9 is the scattering angle defined such that cos 9 = 
r ■ k/(rk) and f(9) is the scattering amplitude. As we 
have shown in Sec. HH this wave function satisfies the 
singular boundary condition at the origin 

ij}(r,k) = ^- [1/r + nlogr] + (3 + 0(r logr) (73) 

47T 

with a I p — —ig. It is also the solution of the integral 
equation (IT4l with A = ig, and consequently due to (fl~5l) . 
(ITBl and (|28p is can be written as 



ip(r,k) = ipc(r,k)-ig 



Y 2 (l + ir ] )e-i"-il 2 W_ lll . h {-2ikr) 



1 + igC(k) 



(74) 

By examining the asymptotic behavior of this wave func- 
tion as T — y oo we can determine the scattering ampli- 
tude f(9). The asymptotic expansion can be obtained 
from (|74[) by using the well known asymptotic represen- 
tations for the Coulomb wave function ipc( r , k) and the 
Whittaker function. The first term in ([74| has the form 
of (f72j) with the Coulomb scattering amplitude [15] 



He) = - 



v 



2fcsin 2 9/2 



cxp{— i?ylog(sin 2 9/2) + 2ictq}. 



The Coulomb phase shift here is given by Co = argT(l + 
irj) . The asymptotics of the Whittaker function as r -> 
oo reads [16| 

; (-2ikr) = e -^/2 e ikr-ir,lo S (2kr) + Q{l/r). (76) 



w 



Introducing these asymptotic representations into ([74]) 
we recover the asymptotic expansion (|72j) and as the re- 
sult we obtain the exact representation for the scattering 
amplitude f(9) 



f(0) = fc(0)+f'(9), 



where 



/'(*) 



ig r 2 (l + irfie-xy 
4^ l + igC(k) 



(77) 



(78) 



The term f'(9) is the consequence of the annihilation 
potential and vanishes if g — > 0. Note that actually /' 
is angular independent, in accordance with the fact that 
the zero range potential contributes only in the s partial 
wave. 

As in the case of the positronium spectrum, this scat- 
tering amplitude can be compared with the result of the 
distorted wave Born approximation (DWBA). The first 
order DWBA correction to fc(9) is given by [l5[ 



f B = - 



1 

4tt 



dr ip c (r, -k)igW(r)ip c (r, k), 



(79) 



where ipc( r ,k) is the scattering Coulomb wave function 
given by (fTTj) . The action of the pseudo potential W(r) 
on the Coulomb wave function can easily be evaluated 

W(r)ip c (r, k) = ip c (0, k)S(r) = T(l + irj) e - mi/2 S(r). 

(80) 
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Thus, the integration in (|79l) results 

f B = _W_ T 2 il + ir))e -,r, 



This, in fact, coincides with the leading term of our result 
of Eq. (|78"1) as q -> 0. 



VII. ANNIHILATION CROSS SECTION 

The total annihilation cross section accounts for the 
loss of flux due to annihilation. This loss is due to the 
imaginary absorbing potential. In the case when H =- 
A r + V\{r) + igVi{r), where V\{r) and V^r) are short 
range potentials, the scattering wave function has the 
asymptotic form 



(r,k) ~ expfir ■ k] + A{6) 



exp[ikr 



The optical theorem connects the imaginary part of the 
forward scattering amplitude to the total cross section 
[lEj]. If an absorption potential is present, the standard 
form should be modified as 14] 



4tt 



3mA(0) - cr 



-9 



(82) 



(83) 



The right-hand side term of (|82|) is conventionally inter- 
preted as the absorption cross section 



k w k 
where a is the total cross section 

a = 2w [ d6 sm9\A(6)f 



(- 



drV 2 (rMr)\< 



(84) 



In case of long range potentials, when Vx(r) is the 
Coulomb potential, the left hand side of the Eq. (|82p 
is not well defined. The scattering amplitude in the for- 
ward direction is infinite and the integral (1831) for the 
total cross section also diverges. The right hand side 
should also be redefined if the pseudo potential is used 
for V2. In order to obtain a generalization of the opti- 
cal theorem (j52")l for Coulomb Hamiltonian with the zero 
range potential we consider the partial wave expansion 
interns of Legendre polynomials Pg(cos6) 



/(0)=]T(2^+l)/^(cos0), 



£=0 



where 



ft = ft + ft- 
Here the Coulomb partial amplitude reads 

f C e 1 

Jt - 



2ik 



(85) 



(86) 



(87) 



with a i = argT(^ + 1 + iri). The additional amplitude f e , 
which is due to the annihilation potential, vanishes for 
I 7^ 0, and for I = reads 



ig T 2 (l +ir))e- n 
"47T l + igC(k) 



(88) 



The problems with the left hand side of the optical the- 
orem (|82p in the Coulomb case come from the strong 
singularity of the Coulomb amplitude (1751) in the for- 
ward scattering direction 9 = 0. This latter singularity 
is the consequence of the slow convergence of the series 
(|85|) pjjj . In order to address this problem let us truncate 
partial-wave expansions 



f N {9)=J2(2£ + l)f i P i (cost 



(89) 



and 



1=0 



N 



a N = ^4tt(2^+1)|/ £ 



(90) 



with N finite nonnegative integer. Introducing the par- 
tial s-matrix elements 



S e = 1 + 2ikf t 



we write the left-hand side of Eq. 



N 



4tt 



(91) 



$mf N (0)-a N .{92) 



If the absorption is absent the scattering is unitary, i.e. 
\S t \ = 1, and (52]) yields 



47T 

— Sm/jv(0) - cr N = 0. 
k 



(93) 



As long as N is finite, this formula is valid for both short 
range and long range cases. The difference arises when 
one is attempting to calculate the limit N — > 00. Indeed, 
in the short range case, the individual limits of each term 
in the left hand side of (|93|) do exist and that leads to 
the standard optical theorem. However, in the long range 
case, the individual limits do not exist, whereas the col- 
lective limit from (|93[) is immediate, and the optical the- 
orem takes the form 

lim ( ^3mjW(0) - a N ) = 0. (94) 

Let us consider now the case when besides a Coulomb 
potential an annihilation potential is also present. Since 
the annihilation effects only the s wave, = 1 for£ > 1, 
and I So I 7^ 1- More precisely, from Eqs. (|86H88[) we have 



So 



e 2lao +2ikf , 
e 2w7f , I > 1 . 



(95) 
(96) 
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This leads for the left hand side of (|92 



JV 



^5>+l)(l-|^| 2 H^(l-|S | 2 ). (97) 



e=o 



With this equation we obtained a variant of the optical 
theorem with finite N 



4tt 



— Qta/^O) - <j n = ^ (1 - |5 | 



(98) 



However, the right hand side does not depend on N, does 
the N — > oo limit exists and yields 



'4-7T, 



JV- 



lim — Sm/jvCO) - ^ = _ (1 - |S | 2 ) . (99) 



If the value of j3 is calculated from this formula pertur- 
batively, i.e. by using the leading term of ip when g — > 0, 
we get 



i>(r,k) = i! e (r,k) + 0(g), 



and 



P = fc(0,k) + O{p). 
Then for the annihilation cross section we have 



^n~y J dr5(r)\yj c (r,k)\ 2 
and for the optical theorem we get 



(105) 

(106) 
(107) 



Um^ ^Sm/iv(0) - ffN \ ~-yJ dr«5(r)|^c(r,fe)| 2 . 

(108) 



As in the case of standard optical theorem (|8"2"j) . the quan- 
tity on the right hand side is the measure of the non uni- 
tarity of the S matrix and in our case it determines the 
annihilation cross section 

<Tan„ = £ (1 - |So| 2 ) ■ (100) 

In order to finalize the generalization (|^|) of the optical 
theorem we need to represent the right hand side of (|100[) 
in terms of the wave function. This can be accomplished 
by using the exact representation (|88|) for the amplitude 
/q and the properties of the Coulomb amplitudes. After 
some calculations, which are given in the Appendix, we 
arrive at the desired representation 

l-|So| 2 = ^|/3| 2 , (101) 

7T 

which leads to the final form of the annihilation cross 
section 



VIII. CONCLUSION 

In this paper we have constructed the zero-range po- 
tential for the two particle system with the Coulomb in- 
teraction. The Coulomb potential modifies the asymp- 
totic behavior of the wave function at small distances. 
This asymptote has been studied with the help of the 
Lippmann-Schwinger integral equation. The pseudo po- 
tential has been constructed such that the wave function 
obeys the singular asymptotic boundary conditions. Al- 
though, the pseudo potential is constructed for the two 
body case, its form is universal, and it can be used for 
several particle systems [201 ]. This work places the delta 
functional model of the annihilation potential | on a 
sound mathematical basis. The optical theorem, which 
formally cannot be formulated for charged particles, is 
reformulated in such a way it becomes applicable to the 
case of Coulomb plus imaginary absorbing potential of 
the zero range. This result facilitates the use of the anni- 
hilation potential for calculating annihilation in collision 
of positrons with atoms, includin g re arrangement pro- 
cesses with positronium formation [14j . 



dann = -y\P\ 2 (102) 

and to the final form of the optical theorem 

lim (^ m f N (0) - a N ) = ^|/3| 2 . (103) 
n^oc \ k J k 

In previous sections the parameter (3 is represented 
through the wave function by the integral 

(3= I drW(r)i/j(r,k). (104) 
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Appendix: Derivation of the representation for the 
annihilation cross section 

Here we present the derivation of Eq. (I101[) . Using the 
exact representations for the Coulomb phase shift and for 
/o w e get 
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1 ' |5 ° l = ^TTi^p 2 7 r|l + z,gC(fc)| 2 (-2S m C(fc) + -|r(l + z 7? )| e 



(A.l) 



I 

We show that the second term on the right hand side 
is zero. Indeed, the first term in figure brackets can be 
transformed as 



2SmC(fc) = -- ^-3m [log(-2ifc) + + ir))] = -J- + ^ - -£-3m^(l + ir?). 

2f Zf Z7T 4 Z7T 



To proceed further, we use the formula [l 

3m^(l + iy) 



+ J coth(Try) (A.3) 



to get 

fc n n 
-2^ + 4"2^ W(1 + i??): 

and therefore end up with 

- 23mC(fc) = 

Now, if we use in the term 
k 



2 V e 27 "> - 1 



n 1 

~2e 2m 'i - 1 



-|r(l + iry)| 2 e-^ 



(A.4) 



(A.5) 



I 

the relation [l6[ 

T(l + iy)T(l - ly) 

we get 



Try 



sinh(7ry) 



A|r(i + i7? )| 2 e — 1 = ~ - — . 

2tt' V " l 2e 2 ^-l 



(A.2) 



(A.7) 



(A. 



So, we have shown that the expression in figure brackets 
in (|A.1[) equals to zero. Thus, we arrive at the equality 

gk\T(l + ir,)\ 2 e-^ _-gk a 

1 ~ 1501 " ^TTWF " ' (A ' 9) 



(A. 6) which proves f|lQlj> . 
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